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D ■ The renormalized energy-momentum tensor of vacuum has been deeply explored many 

years ago. The main result of these studies was that such a tensor should satisfy the conser- 
vation laws which reflects the covariance of the theory in the presence of loop corrections. 
In view of this general result we address two important questions, namely how to implement 
the momentum cut-off in a covariant way and whether this general result holds in the theory 

p; 

with Spontaneous Symmetry Breaking. In the last case some new interesting details arise 



and although the calculations are more involved we show that the final result satisfies the 
conservation laws. 
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I. INTRODUCTION 

Traditionally, the calculation of quantum corrections to the stress-tensor (also called Energy 
and Momentum Tensor - denoted as EMT in what follows) of vacuum is one of the most important 
issues of Quantum Field Theory in curved space-time. The reasons for the special interest is 
this problem are becoming obvious if we remember that the matter fields and particles enter the 
cosmological and most other gravitational equations in the form of EMT of matter, which is usually 
taken as a fluid. The quantum effects of field fluctuations turn out to give some corrections to 
the corresponding equations of state. The most relevant example, probably, is that the EMT for 
radiation gains a non-zero trace due to the conformal (trace) anomaly, which really changes the 
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equation of state for the radiation with some possible relevant effects for the radiation-dominated 
Universe Opposite to this case, the equation of state for the massive particles and baryonic 
matter in general, does not change essentially, because quantum corrections can not make such a 
matter content to be relativistic. 

The situation is quite different in the case of vacuum quantum effects, which can be much 
more relevant than those of the matter sector. Recently there were many publications on this 
subject, including the ones where the possible quantum effects of quantum massive matter fields 
on cosmology and astrophysics were explored. In particular, it was noticed in 

0-Q] that such 

quantum corrections can be defined up to a single free parameter v on the background of general 
covariance. Some observational consequences of the possible quantum corrections were explored in 
0,3] and led to establishing an upper bound on the magnitude of v. Furthermore, the same unique 
form of quantum corrections was applied also to astrophysics |5[ and was shown [l0[] to provide an 
accurate description of the rotation curves for some sample set of disk galaxies, without introducing 
a large amount of Dark Matter content (see for more examples). Some other applications to 



cosmology and astrophysics were also discussed in Ref. 

The applications mentioned above are based on a single, however nontrivial assumption of the 
existence of relevant quantum corrections in the low-energy vacuum sector. Needless to say that 
the most desirable development would be to derive such quantum effects on the regular basis in 
the framework of some rigorous QFT approach. The problem was discussed in [3] and the final 
conclusion concerning existing regular methods was essentially negative. The required quantum 
correction to the effective action of vacuum should be given by a sum of infinite products of the 
curvature tensor components with an infinite number of non-local insertions, hence there are small 
chances for a practical realization of such a calculus. One can note that the situation becomes much 
more definite if we give up the covariance and use, for instance, the conformal parametrization of the 
background metric. In this case it is possible to calculate the quantum corrections However, 
this method is not really safe and is anyway applicable only at the high-energy regime when the 
minimal subtraction procedure is supposed to be reliable. 

It would be very nice to have some alternative approach to the derivation of desirable quantum 
corrections. Recently there were some publications where the result was obtained by means of the 
cut-off regularization in the conformally flat cosmological metric case [if]] (early version) and [n| 
(see also [171]). The main idea is to perform calculations of the "energy density" and "pressure" 
of the vacuum in the momentum cut-off regularization, taking the expansion of the Universe into 
account perturbatively, order by order in the Hubble parameter H . The zero-order approximation 



has been considered before by Akhmedov in [ijj and earlier by DeWitt [uj ]. The output of the 
non-covariant procedure is not the naively expected equation "equation of state" p vac = —p vac 
of the cosmological constant, but the one for the radiation p vac = p vac /3, which led to several 
attempts to understand this result and even to correct it at the ad hoc basis [20I ]. In fact, DeWitt 
explained the result in a very general terms as being produced by the non-covariant regularization. 
The calculations of [15I . [lfj were based on the subtraction of the flat-space result of [la ] , which led 
to the new "equation of state" for the vacuum, this time proportional to H 2 times the square of the 
cut-off parameter. The main problem with this result is that it apparently contradicts either the 
general covariance of the effective action, or the locality of the requested counterterms. in this case 
we meet a violation of the well established fundamental features of renormalization in curved-spaces 
(see, e.g., books 21, 22] and recent papers [3]). However, the results of these calculations should 
be considered as a motivation for the study of possible existence of the 0(-£f 2 )-type corrections to 
the vacuum energy in cosmology. At the same time it looks very important to better understand 
these results at the technical level. This consideration is one of the motivations for the present 
paper. Furthermore, it is interesting to see how the calculations in the cut-off regularization can 
be done covariant. This problem has been recently solved in [2j] on the basis of local momentum 
representation in Riemann normal coordinates (alternatively, one can achieve the covariance of 
finite expressions by imposing the conservation law step by step when adding specially adjusted 



non-covariant counterterms 
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251] ). Furthermore, there is one more possibility which deserves to 



je checked in full details. The cosmological constant term consists of the two main contributions 



261 ]. namely the vacuum classical term and the induced term. The no-go statement of [13l ] concerns 
only the quantum contribution to the vacuum part and, therefore, there is a chance to meet 0(H 2 )- 
type quantum corrections to the vacuum energy from induced part. As one can see in what follows, 
for the induced contribution the route from effective action to the EMT is not so direct as it is 
for the vacuum counterpart. The corresponding calculation requires more efforts and concerns the 
main purpose of the present paper. We shall derive the quantum contribution to EMT from the 
induced term in the covariant way, in the linear in curvature approximation and will eventually 
show that in this approximation EMT of vacuum is local, satisfies the conservation law and hence 
it is given by a linear combination of the metric and Einstein tensor. 

The paper is organized as follows. In Sect. 2 we present a brief summary of renormalization in 
curved space-time and discuss the non-covariant results obtained on the cosmological background 



from this perspective. In Sect. 3 we consider, following 



27] , the spontaneous symmetry breaking 



in curved space and derive the corresponding classical vacuum EMT in the linear in curvature 
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approximation. Sect. 4 is devoted to the conservation law for the EMT of the vacuum in the 
theories with SSB. In Sect. 5 we present as additional technical discussion of the classical EMT 
of the vacuum and its physical relevance in different theories. In Sect. 6 we derive the one-loop 
quantum correction to this EMT. Finally, in Sect. 7 we draw our conclusions and present some 
additional discussions. Some calculations concerning the normal coordinates and local momentum 
representation are addressed in Appendix A and a detailed derivation of equations of motion in 
the linear in curvature approximation is contained in Appendix B. 



II. BRIEF SUMMARY OF RENORMALIZATION IN CURVED SPACE 



The renormalization of quantum 



many investigations starting from 



theory of matter fields in curved space-time was subject of 



281 ] . The most simple way to remove divergences by the consis 



tent renormalization procedure is related to the effective action method 
means of Batalin-Vilkovisky formalism 



22 



29 



301 ] (including by 



231]). The result of all these studies can be formulated in a 
simple form as follows: the theory of quantum matter fields which is renormalizable in flat space 
can be formulated as renormalizable in curved space if there is a regularization which is consistent 
with general covariance from one side and the gauge symmetries of the theory from another one. 
The renormalizability means that the divergences of effective action (at any loop order) are local 
and general covariant expressions compatible with the given gauge symmetries. 

From the effective action perspective the renormalization of EMT is looking quite trivial: one 
has to derive the effective action T and take the variational derivative 

2 (5r 

(Tp,(x)) = g^ a (x) g vf3 {x) —- . (1) 

y/-g(x) og aj3 {x) 

After that one has to introduce the counterterms into the effective action and add them to the 

r in ([1]), which equivalent to performing some very special subtraction of the divergent terms. 

This subtraction should exactly correspond to the covariant and local counterterms in the effective 

action. After that the coefficients of the remaining finite terms should be fixed by imposing the 

renormalization conditions on the renormalized classical action and/or renormalized EMT. For this 

end such a classical action should be chosen in a special way and include all the structures which 

are possible to emerge as counterterms. 

The arguments based on covariance, locality and power counting lead to the following form of 

the classical action of external metric (vacuum) : 

Svac = Seh + Shd , (2) 
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where Seh is the Einstein-Hilbert action with the cosmological constant 

Seh = J d 4 x^(R + 2A) . (3) 

and 

S H d = J VaV=#{aiC 2 + a 2 E + a 3 OR + a 4 R 2 } . (4) 

Here C 2 = R* vap - 2R 2 j3 + (1/3) R 2 is the square of the Weyl tensor and E = t& p - 4R 2 a/3 + R 2 
is the integrand of the Gauss-Bonnet topological term (Euler density in d = 4). Let us remark 
that the presence of higher derivative terms and cosmological constant are necessary to have a 
renormalizable theory. 

In the present paper we will be interested to perform covariant calculations around the flat 
space-time in the linear in curvature approximation. This means we will systematically ignore the 
higher derivative part @ and, in general, will not pay attention to the 0(R 2 ) and 0(ni?)-terms. 
This means, in particular, that the form of the divergent structures which one can meet in {T^ u ) 
is restricted to the two terms, namely the ones proportional to g^ u which are responsible for the 
renormalization of the cosmological constant term and the ones proportional to the Einstein tensor 
G^ u = R^u — (l/2)Rg flu and responsible for the renormalization of the Einstein-Hilbert term in 
the effective action. 

Our calculations will be performed in the local momentum representation, based on the use of 
Riemann normal coordinates. Also, we shall use very simple cut-off regularization in the Euclidean 
local momentum space. This regularization has been shown equivalent to the cut-off of the pro per 
time integral in the Schwinger formalism in flat space 



311 ] and recently has been used in 



32] to 



calculate effective potential of the scalar field in curved space-time. 

An alternative approach to renormalize EMT in curved space-time is to work directly with the 
classical expression for the EMT and perform calculation. This approach is the most traditional 



one sec 



21 



1 and references therein). The covariant calculations in this way have been performed 
34J by means of the point-splitting regularization, without or with the use of effective 



in [33J and 

action method. The covariant structure of divergences of EMT which has been described before is 
restored in the limit of zero split ting , but only if this limit is taken in a special invariant way. 
Let us consider the result of 



33| for the quartic divergent part of the quantum corrections to 
EMT. For the sake of simplicity we can deal with the flat space expressions, because the quartic 
divergences are not really affected by this choice. Then 



11/ . n^ny 



< T ^>quart. div " 2tt 2 n a n« \ 9fiu 4 nprtP) ' (5) 
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where n a is a small non-null four-vector denning the point splitting regularization of the corre- 
sponding Green functions G(x,x) — > G(x,x + n). Now, if we chose the vector n in temporal 
direction n = (e 2 , 0.0.0) we get that 

/3 0^ 



1 



quart, div 



2vr 2 e 4 



10 
10 

y o o o i ) 



(6) 



which is a traceless quartic divergent component of the total energy- momentum tensor. This 

result is in accordance with the one based on naive momentum cut-off in flat space of [18fl. As it 
□ 

was explained in [19|], there is no contradiction with the expected local Lorentz invariance of the 
divergences, because the origin of the © is in the use of a non-covariant regularization. In case of 
the point-splitting with temporal direction the breaking of Lorentz invariance is due to the non- 
relativistic choice n = (e 2 , 0.0.0). In case of cut-off regularization the origin of a non-covariance is 
different but since it is equally non-relativistic, the final result is the same. 
One expects that a Lorentz invariant regularizations would give rise to j3] 

/ 1 \ 



1 



v^iw) quart, div 27r 2 e 4 



(7) 



0-10 
0-10 
\ -1 J 

which is proportional to the Minkowski metric tensor and can be interpreted as a standard divergent 
contribution to the zero point energy or cosmological constant term. 

Let us show how this occurs in Pauli-Villars regularization [35] . The contribution of a free scalar 



field with mass m to the vacuum energy p is given by 



i 



16vr 2 



8 2 m \\l J 



(8) 



where Q is a 3-momentum space cut-off. Whereas, the same contribution to the pressure reads, 

1 



P 



- m 2 n 2 



-m 4 + -m 4 log — + O I — I 

8 2 m 



(9) 



48vr 2 V" "' " 8"' ' 2"" "° m ' \S~lJ 
The Pauli-Villars regularization is defined by a family of scalar and ghost fields with masses 
/x?M 2 + m 2 , i = 1,2, - ■■ N with degeneracies Sj. Positive degeneracies correspond to scalar 



fields and negative degeneracies to ghost fields. The Pauli-Villars conditions 

N N N 







(10) 



8=1 



i=l 



i=l 
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guarantee that for a free field theory with mass m the quantum corrections to the vacuum energy 
and pressure are finite, i.e. all f2 quartic, quadratic and logarithmic divergences are canceled out. 
Notice that there are always non-trivial solutions of the Pauli-Villars conditions equations (110p . e.g. 
s = (1,1,-2,-1), /j 2 = (5,8,2,9). However, in the limit when the mass M of the Pauli-Villars 
regulators goes to infinity we recover the quartic divergences of the vacuum energy-momentum 
tensor which now are of the form 

/ 1 \ 

cM 4 



M"/ quart, div 



mV 4) 



2tt 2 



0-10 
0-10 
\ -1 J 



(11) 



where 



1 N 



i=l 



is an arbitrary constant given in terms of the regulating Pauli-Villars parameters. Notice that the 
sign of the vacuum energy correction might become positive or negative depending of the choice of 
the regularization. In spite of the use of a non-covariant auxiliary cut-off the final result is covariant 



361 ] . However, the existence of an ambiguity in the leading quartic divergence and its sign is a 



puzzling characteristic of quantum vacuum. From a renormalization viewpoint the ambiguities 
can be traced back to the locality of the cosmological constant term of the effective action. The 
effective value of the coupling has to be fixed by an explicit choice of renormalization prescription. 

In the case of sub-leading divergences something similar occurs. The quadratic divergences also 
acquire a covariant form in Pauli-Villars regularization 

^10 o ^ 

c'M 2 m 2 



V^Wquad. div 



M 2 t$ 



2vr 2 



0-10 
0-10 
\ -1 J 



(12) 



with 



1 N 



(13) 



i=l 



In the case of quadratic divergences there is one more specific ambiguity, especially when they are 
calculated in the special cosmological background depending on the Hubble parameter H. Imagine 
we have obtained the result in the form 



<T M „) = MH$ + M 2 H 2 t$+ .... 



(IV 



(14) 
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Now, in this expression H is effectively used as a constant, and therefore we can redefine the cut-off 
as M 2 —7- M' 2 = M 2 + XH 2 , where A is an arbitrary dimensionless parameter. As a result we 
arrive at the new form of the power-like divergence, 

(V) = + M 2 H 2 [tg + 2f$] + ... . (15) 

with even greater degree of ambiguity. In the case of a theory of the quantum field with mass m 
one can perform a more general redefinition M 2 -)• M' 2 = M 2 + XH 2 + Tfli , with even more 
ambiguity, etc. It is important that the logarithmic divergences are not affected by this ambiguity 
and, in general, represent the most universal and well-defined part of quantum corrections [37J (see 



also 38] for a recent discussion of the subject and further references). 

One simple way to get free of the mentioned ambiguities is by using the effective action method. 
The prescription which we have already described above is simple. First one has to derive the 
divergent and finite (at the level which is possible) of effective action, add counterterms, perform 
renormalization. At the second stage it is necessary to take a variational derivative with respect 
to the metric (pQ) and obtain the divergent part and/or renormalized EMT. In the next section we 
shall see that this procedure works even in the situation with SSB, where the procedure described 
above is essentially more complicated than in the free field case. 

III. SSB AND EMT IN CURVED SPACE 



We start following Ref. [27f] . However, since our purpose is to consider the most simple model 
with spontaneous symmetry breaking (SSB) in curved space, we will consider the single real scalar 
field, while in the mentioned reference the charged scalar was used. The classical action of the field 
(j) with a non-minimal coupling and a self-interaction is 

S sc = I d A x^-g^-g^d^dA + \m 2 ct? + ^R<p 2 - . (16) 

The dynamical equation for (f> has the form 

- + m 2 + ^ - -A0 3 = 0. (17) 
6 

Consequently, the vacuum expectation value (VEV) for the scalar field is then defined as solution 
of the equation 

-n v + m 2 v + t:Rv-lxv 3 = 0. (18) 
6 
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It is easy to see that there is no constant solution for this equation for £ 7^ 0, while the value £ = 
is inconsistent with renormalizability of the theory (see, e.g., [22]). Hence we can find the solution 
for the vacuum expectation value v only in the form of the power series in £R, 

v(x) = v + vi(x) + v 2 (x) + ... . (19) 

In the zero-order approximation we meet the conventional flat-space expression, 

4 = ?£. (20) 

As we have already mentioned above, in this paper we will use the approximation of small curvature 
and are interested in the first-order approximation only. At this level one can easily find a non-local 
expression 

Vl = (21) 

In a similar way, it is possible to construct further approximations, but this is beyond the scope 
of the present paper. Thus, let us concentrate on the expression (|2ip and simplify it further by 
neglecting terms with derivatives of the scalar curvature. Of course, this approximation works only 
for the sufficiently large value of the square of the physical mass of the scalar excitation near the 
point of the minimum, 2m 2 = Avq/3. Then we arrive at the quantity 

Vl « P-R. (22) 



It is clear that the same solution can be obtained directly from Eq. (|18p if we disregard the term 
Dv and use (|20p . In our opinion the approach followed here is better, because it enables one to 
control the approximation. In further calculations we shall use the expression for classical solution 
of the theory in the point of the minima of the SSB problem, 

00c = v = vo + vi , where Vq = ^™ and v\ = p- R . (23) 

The renormalization of the vacuum sector of the theory with SSB has been described in great 
detail in [27!], thus, we shall not elaborate on it here. Instead, let us discuss the definition of the 
EMT at quantum level for the theory with SSB. The classical energy-momentum tensor of the field 
<fi in the external metric field g = is defined by the relation 

2 5S[g,4>] 

= — g m g v p — . (24) 

V~9 °3a(3 
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At quantum level, the energy-momentum tensor is given by 



hi/ i 



-g 



QliaQvP { 



6S[gJ] 



(25) 



where <j> is quantized field, <fi ~ ua) + u*a and a |0) = 0. As far as g^ a is classical external field, 
so we can take it out of the ( | . . | ) freely. 

According to our previous discussion, we will follow the functional representation of Quantum 
Field Theory, where the basic object is the generating functional of vertex function, or effective 
action, T = T[g, For the case of a scalar field it is defined as a solution of the functional equation 
(see, e.g., 22| for a introduction) 

sr\g,<f>] 



ex P 1 fr T i9,< 



exp 



+ 



(26) 



In this work we restrict consideration by the one-loop approximation, when the effective action 
in Eq. (pQ) becomes the sum of the classical term and of the one-loop correction 



r<%, = S[4>, g^} + «T«[0, g^] . 



(27) 



Then the one-loop EMT of the vacuum can be cast into the form 



(T^(x))« = T^(s)+T«(s), 



(28) 



where the first term is classical contribution, 

2 



-g 



g^ia gu/3 



5S 



and the second one is one- loop correction to it, 

2h 



■ gfia guf: 



8g a p 
«5f« 



<i>-Ht>o 



3g a /3 



<j>->4>o 



(29) 



(30) 



In both cases 4>q is the solution of the equations of motion. If we deal with purely classical theory, 
then one has to replace in (j29|) the value <po = <fro c from Eq. (|23p . After the one- loop correction is 
taken into account, we have 



5S[g,M | h 6jS% 1 0o\ = Q 



(31) 



where the replacement <fi — > (j)o should be performed after variational derivative. The Eq. (I3ip can 
be solved by iterations in h. At one-loop level 



+ Hi > 



(32) 
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where 4>q c is the classical solution f|23|) . In order to find <j>\ one has to replace f|32|) into (|3ip . In 
the first order in h we meet the equation 

* 2 %.« #1 + £F%M =0 , (33) 



and obtain the solution in the form 



(34) 



and, therefore, 

00 = <foc - M x!TZ ULJ ) ~~ — ET yULJ • ( 35 ) 



One has to replace this formula into the expression for EMT, 

{Tnu) = j= 9 m g v p y^— \s[g,4> ] + ht (1) [c/,</>o]) • (36) 

In this way we arrive at the general expression for the EMT in the scalar theory with SSB, 
,rr \ 2 )8S[g,<f>0c] 

{J-fiu) — -j=giia9vP 



-9 { og a/ 3 

+ j. STWfafa] _ h 5 2 S[g, Oc ] / 5 2 S[g, Oc ] \ ~* 8 fW [g, Oc ] 1 

The first term inside the brackets is classical, the second is typical for the free theory and actually 
does not depend too much on the kind of such theory. The last term emerges only due to the fact 
that we deal with the interacting theory. In the free theory this term is zero. 
It proves useful to define 

(T^) = (T M „}„ + (T^i, (38) 



where 



and 



{T^) v = -=g m g v $\ — - Vn > (39) 

V^g 5g a p $g a /3 



lrr v 2^ j s 2 S[g,M [ 8 2 S[gAoc} \ ' $ [g , M ■ , m . 

{l^vji = -j=g m guii < ^- ^ I — — I — f • ( Kb 



The quantities (T^ u ) v and (T^ v )i represent the vacuum and induced parts of the EMT, respec- 
tively. Both quantities will be calculated in this work in the context of SSB. 
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IV. COVARIANCE AND CONSERVATION OF VACUUM EMT 



The conservation of the energy-momentum tensor of vacuum is always regarded as the main 



requirement for the consistency of the theory (see, e.g., 2lj and further references therein). As far 
as we deal with the lower-derivative approximation, the condition of conservation V^(T flu ) = 0, 
together with the requirement that the EMT should be derived as variational derivative of covariant 
effective action, can fix the algebraic form of {T^ u ) completely, leaving the room for only two 
numerical parameters in case of quadratic and logarithmic divergences and finite part and for a 
single numerical parameter for the quartic divergences case. 

The reasons for this special importance of the conservation law are as follows. For the divergent 
jarts of effective action the situation is especially simple, because we know it should be local (see 



231 ] for a recent discussion of this issue in curved space-time). As we have already explained in the 
previous section, in the lower-derivative sector this means that the possible counterterms have the 
form of the Einstein-Hilbert term and of the cosmological constant term ((3|). Consequently, the 
divergent part of the vacuum EMT should consist of only two structures, namely 

(T,u,) = C x g^ + C 2 GV, (41) 

where C\ = A^fi 4 + /c2^ 2 + &l In (^/Mo) an d C*2 = + II In (^V^o); with ^4, /C2, &l and h, II 
being numerical constants. The values of these constants depend on the choice of the quantum 
theory, on the order of loop expansion, but the structure of divergent part must be always like in 

Concerning the finite part of the EMT of the vacuum, it is possible to have much more com- 
plicated expression than the one presented in (|4ip . as a result of resummation of the series in 



curvatures and Green functions 
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One can have an indication to this possibility, e.g., from the 



calculation in conformal variables [141 ] . However, as far as we are going to perform a relatively 
simple calculation in the 0(i?)-approximation, there is no room for non-localities in the effective 
action, so what one should expect as a result is the same expression (|4ip . 

Needless to say that (14ip is the only form of EMT which can be derived from some action 
principle and also is the only form which satisfies conservation law. Let us start from a well-known 
derivation of this relation in general case and then consider the same thing in view of Eq. (|37|) . 

The effective action T is covariant scalar functional depending on metric g^ v and scalar field (j>. 
If we perform infinitesimal general coordinate transformation x a —> x' a = x a + £ a (x), these two 



13 



fields transform according to the known rules 

Sg^ = -V„£„ - , 8<j> = - ed^4> . (42) 
Then the identity corresponding to diffeomorphism invariance of V is 



Now we take into account that the functional derivative vanish on-shell, that means 

^%M = . (44) 

Then integrating the first term in (|43p and taking into account the definition of EMT ([1]), we arrive 
at the conservation law, V M (T M1/ )0 o = 0. 

Let us now see how the same considerations look when we perform the expansion of T[g, cf>o] 
into series in h. Here we are interested in the expansion up to the first order and the main question 
is whether the two parts EMT, namely vacuum and induced ones, (|40p and (|39p . do satisfy the 
conservation law separately or only when they are summed up. 

At zero order everything is quite obvious, for we have 

M = , and V^ L =0. (45) 

At the first order in h we notice that the conservation law is satisfied only on-shell. Now, the 
solution cj)Q = (fto c + fufti of (|34p , was found exactly to provide that 4>q is the solution of the effective 
equations of motion at one loop. Hence, we should expect that neither one of the two terms (|40p 
abd (I39p will satisfy the conservation equation and only for their sum this equation must be valid, 

V^T^oc)), + V^(V(0 Oc ))i = 0. (46) 



On the other hand, this means that the sum (|38p should have the form (I4ip while each term can 
have more arbitrary form, for example the Ricci tensor and scalar curvature term may not form 
Einstein tensor. 

Finally, both the covariance arguments and conservation law indicate that the quantum EMT of 
vacuum, in the low-energy sector of the theory, must have the form (j41h even in the presence of SSB 



which produce much more sophisticated forms of EMT, such as (I37D . The restricted form of the 
vacuum EMT (|4ip should hold even at higher loop orders, at least for divergent contributions [42]. 
This is a strong statement and it is worthwhile to check it by direct calculation, at least in the 
one-loop order. We shall do it in the next section. 
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V. DERIVATION OF STRESS TENSOR: CLASSICAL PART 

In this section we shall derive the EMT of vacuum at classical level and then, in the next 
sections, turn to the one- loop contributions. 

Let us perform the calculation of the EMT T^ u in the vacuum state, which is characterized by 
the VEV of scalar field defined in Eq. (|23p . The calculation of T^ v is rather trivial and we obtain 

T, v = --L^^^M = (2^-i)^(V0) 2 + (l-2O(^)(a,0) 

+ 2^(^,00 - V^V^) + £<£ 2 (iV - + -^mV + Y± 9 ^ ■ (47) 



The trace of the scalar EMT on-shell (|17p can be easily reduced to the form 

T£ = m - 1) [(V0) 2 + R<P - \ 4 ] + 2(1 - 30 m 2 ^ 2 . (48) 

We observe that for m 2 = and £ = 1/6 we have Tff = 0. However, we are interested in the massive 
case given by Eq. (|23j) . At this point it is worthwhile to discuss the practical realization of the 
O(R) approximation, which we will follow in this section. The main question what to do with the 
derivatives of 4>o c = Vq + v\. Since vq is a constant, its derivative is obviously zero. Furthermore, 
a derivative of v\ gives us 

V a foc = = T- 1 - VaR (49) 

Xv 

and, therefore, goes beyond the limits of our approximation. As a result we can always treat R 
and v\ as constants, that leads to great simplification of all calculation. 

Replacing (|23]1 into (|47p and keeping only terms linear in curvature tensors, after small algebra 
we arrive at 

1 Xv^ 

T^oc) = tv 2 (iV - - Rg^ - 9iiv . (50) 

This expression is nothing else but the induced contribution to the Einstein equations. It is natural 
to attribute it to the gravitational part of these equations, which can be, eventually, written as 
1 1 \ /„ 1 



+ g^— ) (V - - Rg,u) - (pT + P t d ) 9,u = T- atter , (51) 



> 8irG vac 8-rrGind 
where 

_l_=-£«8 and Pt d = -^-- (52) 

In this equation G vac and p v £ c denote the vacuum Newton constant and the cosmological constant 
density, which are independent parameters that are originally present in the action of the theory. 
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Contrary to that, Gi n d and p™ d are induced quantities which depend on the details of the quantum 
theory of matter fields under consideration. 

The induced contributions here are due to the SSB, and an equivalent mechanism of their 
generation is working also for the Standard Model (SM) and Grand Unification Theories (GUT's). 
The values of induced and vacuum cosmological constants are known to be, at least, 55 orders of 
magnitude greater than their sum in (|5ip . that gives rise to the cosmological constant problem 26] 



(see also references therein and [2j in relation to renormalization of the cosmological constant). On 
the contrary the relative magnitude of Gi n d , namely 

Gjnd 8<^q ^ 

^vac 1V1 p 

is small for the SM case where Vq ~ 10 5 GeV 2 . Even if the value of £ corresponds to the Higgs 
inflation, £ ~ 4 x 10 4 , the Planck suppression is strong due to the relatively huge value Mp ~ 
W^GeV 2 and hence the contribution of (|53p is irrelevant. 

However, the situation can be quite different in GUT's, where (in the super symmetric versions) 
we have Vq ~ 10 32 GeV 2 . Then for the mentioned above magnitude of £ we arrive at the estimate 
G{/G v ~ 1 and the effective sum in f)51 1) becomes close to zero. This means the value of 1/G V must 
be taken about twice larger than the observed sum. Hence, the classical screening due to induced 
value may be relevant in this case. In the rest of this paper we shall check that the quantum effects 
do not break the structure of (f5"T]) and calculate quantum corrections to the quantities Gi n d and 
pf d in (52}. 



VI. ONE-LOOP CALCULATION IN THE O (^-APPROXIMATION 

Let us now perform quantum calculations using the expressions for vacuum and induced parts, 
(|40p and (|39p . The calculations will be done in the local momentum representation and covariant 
momentum cut-off regularization. For better organization, this section is divided into subsections. 
First, we consider some general notions, then derive the flat-space result, then present some minimal 
mathematical tools for the local momentum representation, and finally perform derivation of the 
more complicated, curvature-dependent part. 
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A. General considerations and derivation of (7)j„)„ 

Our starting point will be the one-loop effective action, T^[g,(j)]. By construction, this is the 
effective action in the theory with unbroken symmetry. One can write, using derivative expansion, 

T^[g,cf>) = J ^^{-F e// (0) + iv^-^(-^)v^ + ^ 2 %(^)i?+...}, (54) 

where the effective potential part has the form 

V eff (4>) = Vo + ViR + OiR 2 .), (55) 

which was recently calculated using covariant momentum cut-off in [3^] and k^l^^j and k^(^_ 
are the form factors which also contain different powers of derivatives. The expansion in (|54p is 
infinite, but we can easily set the limit on it, following the same approach which was used in the 
previous section. For £ = we know 4>q c = const, according to Eq. (I17j) . Therefore, any derivatives 
of 4>o c are actually proportional to £ and hence to R. 

As far as we are interested only in 0{R)- terms, we can take only constant part of the form- 
factor k^ in (|54p . and also strongly restrict k^ form factor, also by taking its constant part. Hence 
we can trade 

and ^(^)^ X( ^' (56) 
Futhermore, in the given approximation the term |</> 2 x(</>) R is a part of the effective potential 
V\ = V\{4>). So, for us f^ 1 ) [<?,(/>] becomes 

fW(g,<f>) = J d 4 xV^{ -y e// (0) + ^(<A)(V<A) 2 }, (57) 

with V e jf = Vo((j)) + V\((j)) R and </>—>• 4>o c . 
Let us consider 

J tfxJ^gZifflVt) 2 = J d 4 x^-gV^- f d*xJ=gZ(</>)<l>n<j> 

- [ d 4 z,/^'(0)(V0) 2 . (58) 



For (ft — > (f)Q C , the quantity □</> can be written as 

□ + 2m 2 ~" ~ 2m 2 



□0 Oc = U{vo + t*) = — ^ a R = a R + 0(0 2 R) . (59) 



On the other hand, 

V7 rh„ — V7 n„ _l_ V7 iu — V7 nu — 

2m 2 



V M <Aoc = V„«o + V^i = V„«i = |^V M i? + 0(V 3 R) . (60) 
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It is now obvious that, because of □ <C m 2 for c/>o c , the whole quantity (Vc/>) 2 is beyond our 
approximation 0{R). Finally, we can restrict our consideration by the effective potential, using 
the expression 



The renormalized expression of the potential is 



d 4 Xy/=gV e ff((l>Oc) ■ 



22J 



v e r f j(g^, <p) = vr + vr n R 



(61) 



(62) 



1 



1 



1 



2(4«y& V "-^ 2 -(t-i) R ( V "-^ 



In 



V" 



m 



In this expression we used a general form of classical interaction term V = V((p), but later on it 
will be replaced by V = Ay? 4 /4 

For the sake of completeness we will also consider the divergent part of the non-renormalized 
potential, in the local momentum cut-off regularization. In the given approximation we have 
V e f f (g^, <p) = V div + Vf v R, where 



Vr 



div 



div 



(63) 
(64) 



It proves useful to introduce a notation for the one-loop contributions to the equations of motion 
for a scalar field, 



f(1) 



^div ' fin 



1 ST™ 



-9 



dV 



(i) 



(65) 



After adding the corresponding counterterm, we will also have e^en- Let us now remember that 



V eff = V (<f>) + V 1 (<f>)R, 



-r(l) 



vq + vi , and define also 



8Vr 



(1) 



(i) 



dVr 



(1) 



4>oc 



R 



4>oc 



R 



(1) 



--(1) 



+ £ 



(1) 



"0 



(66) 



Obviously, both £q and e$ have finite and divergent parts and after adding counterterms we can 
also define their renormalized versions. In this paper we will calculate divergent and renormalized 
quantities only, but the original finite parts can be calculated in the same way using, e.g., the 
effective potential from 32|. The last relevant observation is that, within the 0{R) approxima- 



tion adopted here we can treat all versions of £q and (divergent, finite non-renormalized, 
counterterms and renormalized) as constants. 
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Let us now calculate the simplest quantum term (T^ u ) v , denned in ([39]) . Starting from (f6T|) we 
can easily arrive at 

\-*-fJ,h>/V — i QllQiQvfi 



3 uu cj u iJ c 

-9 og a/3 

( 1 \ SV*^ 

■2Wi(u )^--i25^J + ftVo(«o)fly«/ + hv i9vv 



= - 2 HViivo) + hV {vo) 9ilv - ^ i?4 1} 5^ • (67) 

This formula is remarkable, because it confirms what we have anticipated in the previous section. 
The first two terms in the last expression are quantum contributions to the Einstein tensor and 
cosmological constant part in the Einstein equations. However, the last term looks odd, for it 
violates covariance, conservation law and can not be derived from the action principle. So, we 
should hope that it will cancel with the corresponding contribution from {T^ u )i in (|40p . which is 
the last term in (|37p . Let us see whether this really happens in the next section. 

B. Calculation of (T^ v )i 

Our first step will be to rewrite the expression (|40p for {T^ u )i in a more useful and detailed form 



{T^uix^i = 2hg fia (x)g u/3 (x) / d A y^J-g{y) / d A z\/-g{z) 



fS\gM\ 



9 ^9ap{x)54>{y) J 

l 

(68) 



/ 1 5 2 S[g,4> 0c ] \ x ( 1 5fM\gM \ 

Let us note that the metric-dependent quantities are always understood through the normal coor- 
dinate expansions (see some details of this technique in Appendix A). 

The next step is to derive all three factors inside the integrals of the Eq. (168p . We need to 
perform this calculation in the 0{R) approximation, which we follow here. The first factor can be 
obtaned by variating ([47]) with respect to <j) or just taking a second variation of the action. After 
some algebra we arrive at 

1 S 2 S 

^(V M V„-(7^n) + (2^-l)(V M 0)V v 



\T I 9 S(/)(y)Sg ftv (x) 

+ (~-2£)^(V^)V A + £(V M V^- 5/iI/ n< 
- &(Rij,v - ^Rgn^j + \^?4>9p,v - q^Oimv ■ (69) 

Now one can replace in the last expression <j) —> (pQ C = v o + v\ and remember that all derivatives 
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of 0oc are beyond our approximation. In this way we obtain 



1 S 2 S 



-g 



$9c 



£4>0c (V M V„ - g^u) 



0C I RfjLU 



4>0c 9iiv ■ 



(70) 



Finally, we replace here the expansion ([T9"]) up to the first order in curvature, 4>q c = vq + vi, with 
and v\ taken from Eqs. (|2Q|) and ([2ip . Also, we use the normal coordinates expansion of the 
operator V^Vj, — g^O which is calculated in the Appendix A. The final result for the first factor 
inside the integral in Eq. (|68|) has the form 



1 



SS 



-9 §9v,vH 



90c 



(71) 



+ 3^0 



2R\ tlu)T y T d x + 2 r/^ R^y T d x + R^y a y p d 2 + v« P «/l/V^ 



The first term in the r.h.s. is of the zero order in curvature and the rest of the terms are of the 
first order in curvature. Indeed, after all calculations are completed, we will trade the metric rjnu 
in the point P to the general one g^, but for a while it is better we write it in the way we did. 
Let us now consider the second factor inside the integral in Eq. (|68|) . 

-1 

= G(y,z;0 Oc ), (72) 



1 SS 2 [gM 



-9 



y.z 



This is nothing else but the propagator of the scalar excitations near the point of the minima. It is 
im por tant to remember that we will need the dependence on the curvature. Therefore, according 
to I39I . I40I (see also [32I ) one has to modify the ([72"]) to the form 

\ -1 



1 



-g(y)] 1/4 l-g(z)Y /4 8<f>( y ) 5<f>(z 



■ : I = G(y,z; 4> 0c ) . (73) 



32, QQ, 



Now we can use the known result for the propagator from the mentioned references 
but first we have to evaluate the mass of the scalar excitations near the point of the minima. 

One can start from the full propagator with <p oc = vq + v\. Starting from the equation ([T7]) we 
arrive at 



1 5 2 S 



-U + m 2 + llR-^4> 2 



(74) 



Next we replace 



4> -)• 4>l c = (v + vi) 2 « Vq + 2v v 1 



(75) 



20 



Replacing (|2ip and (|20p into (|75p . after some small algebra we arrive at 

1 5 2 S 2 f 6m 2 \ 



y/=g 5<j)8(j) V □ + 2m 2 . 

« -(□ + 2m 2 ) - 2£R, (76) 

where at the last step we used the (^-approximation, as it was already discussed above. Now, 
after we compare the last expression with Eq. (117p , it is clear that (176p means we have a propagator 
of a scalar particle with positive mass 2m 2 and with the non- minimal parameter — 2£. By using 
the general expression (|107p we obtain the Euclidean version of the second factor inside the integral 
in Eq. (I68D in the form 



G(z-y) = /^e*(—) 



f U R 



(77) 



k 2 + 2m? V s 6/ [k 2 + 2m 2 ) 2 _ 

The third factor of the integrand in Eq. (|68p is nothing else but the effective equation of motion 
(|65p . According to (|66p we can write it as a sum of classical and quantum parts, e = + he^\ 
where the last can be also expanded into series in scalar curvature, e^ 1 ) = ei + . For the sake 
of completeness we have calculated these expressions, but since they are rather cumbersome, we 
postpone them to Appendix B. 

Now we are in a position to derive ^Tav)i- As a first step we obtain the flat-space expression 
and then consider a bit more complicated curvature-dependent terms. 

In the flat-space limit we have only first terms in the r.h.s. of Eq. (|7ip an d (|77p . and also need 
only £q -parts in the equation of motion (both divergent and renormalized versions). In this way 
we arrive at the expression 

<T^(x)>° = 2ft£«o / d 4 zd 4 y 5\x-y) 



/j4ju p ik(y—z) 

where the upper index indicates flat-space limit. After performing integration over y and using 
(z) = const, we get 

/f rfik P ik{x—z) 
d*z(d„d v -r,»„d 2 ) x j (27r)4fc2 + 2m2 4 1} (^) 

= 2H v f d*k 5\k) *ffi^V e ito = o . (79) 



Thus, the contribution of the last term in (|38p to the induced cosmological constant is zero. 
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As a first byproduct we also obtain that in curved space-time the contributions of the third 
factor, being it &j~l iv or £^il en , are also vanishing. The reason is that both are constants in the 
O(R) approximation and we did not use an explicit form of a constant £q in the calculation 
presented above. 

As a second byproduct we can see that in curved space-time the curvature-dependent contribu- 
tion of the second factor (177p vanish too. The reason is that, if we trade 

1 1\ R 



/ 1 \ R 

( 2e " e) fP + 2m 2 ) 2 (80) 



k 2 + 2m 2 V s 6/ (k 2 + 2m 2 ) 2 

in (|78p . the zero output of the integral will obviously remain the same. So, after all we need to 
take into account only the curvature-dependent terms in the first factor, Eq. (|7ip . 

The last step is to perform the curved - space calculation in the 0(R) order. Taking into 
account the arguments presented above, we arrive at 

/(' f]^\c ^^i.y~ 
dVz / (2^^°^ (y)54( "" y) ^T^' (81) 

i=l 

where 

Oj 3 J = I £v [R x {tlv)T + n^R*] y T dx , 
0$ = ^voR^y^d 2 , 

O® = I Hvo Vlu , R a % y a yPd p d a . (82) 

Let us evaluate all the terms of (|8ip . indicating the term in (|82p by the left upper index. 

The contribution of 0$ has the form which strongly resembles (|78p and can be treated in the 
same way, 

= - 2K vo R,u 4 1} / d*k 5\k) = ~^ R ^ 4 1} • (83) 

(2) 

Next, the contribution of O^J vanish, for it has the same structure as the flat-space term (J75J), 

(2) <T/*>* = 0. (84) 
The contribution of (Oy^J can be presented in the form 

{3) (T/u/)} = ^ [R X w + r,^] 4 X) II , (85) 
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where 



Ml, ft ik(y-z) 



The last integral can be calculated by elementary means to give 



and hence, after some small algebra, we obtain 

(3) <«! = -^(^-^) £ f). (88) 
The contributions of (O)^D and {O)^} can be expressed in the form 

{4 ' 5) (T, u )\ = £ f) [i?^ ^ - i? a ^ CT ^] J p(T , , (89) 



where 



j4jl a2 p ik(y—z) 

' " ' (2vr) 4 v y ' y y dyPdy cr k 2 + 2m 2 v ; 



Taking this integral we obtain 



J pa ^ = ^ <V ^ where 5 pa ^ = \ - $8°) . (91) 
Using this result, after small algebra we arrive at 

(4 ' 5) <V>! = + (92) 

The total expression can be obtained by replacing the results for all contributions ([83]) . ([83]) . 
([55]) and (f92|) into (|8"2"|). After a some algebra we come to the result 

(V), 1 = ^ 4 1} (-R»u + Rn^) ■ (93) 

Obviously, this expression is different from and therefore it violates covariance and conservation 
law. However, if we sum up with the previous result {Tu V )y from Eq. (|67p . we arrive at the 
expression which agrees with our expectations, 

(Tfj,u) = {T pl /)i + 

h ^ V (t> 1 o„ W 1 ) 



2 y i( u o) + |5 4 1} 1 + hV oM 9nu , (94) 
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where we finally replaced the flat metric r]^ u by the general one g^. 

In order to rewrite the quantum contribution in the final form, one needs the expressions for 
^o(^o)) Vi(vo) and e . The renormalized and divergent versions of the first two can be obtained 
from (|62p and (|63p in the form 



V ren (v ) 
V dtv (v ) 



1 



■ m 4 In 



(47r) 2 V ^ 

Tft^ r / Op 1 

jt 3ft 2 - 2m 2 In — = 

32 vr 2 L Vm 2 



(95) 
(96) 



and 



m 



(4vr) 2 
1 



In 



2m" 



32 vr 2 6 



Q 2 + 2m 2 In 



n 2 



m- 



(97) 
(98) 



Taking into account in (|94p the expressions for gg derived in Eqs. ()115[) and (|116p of Appendix 
B, we arrive at the final result for quantum contributions to EMT, 



ft 



fiv /ren 



hm 4 ^ , {2m 2 
In 



(47T) 



m 



(47T) 



2(1 + 30 In 



2m" 



+ H 



G 



[IV 



(99) 



for the renormalized expression and 
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div 



Km? \ 9 9 / SI 2 

3Sl 2 -2m 2 a ^ 

\m 2 



32 7T 2 



,2 



+ i^( 4 «4){ n2 - w "4K < ioo > 



for the divergent one. 



VII. CONCLUSIONS 



We have considered several aspects of the Energy-Momentum Tensor (EMT) of vacuum in 
curved space-time. A naive calculation using a momentum cut-off produces a result which appar- 
ently violates general covariance. It was noticed long ago that this is the effect of the non-covariant 
cut-off scheme and therefore can be hardly regarded to be a physical feature of the theory. The 
two questions naturally arise in this respect, namely whether it is possible to introduce a cut-off 
on a covariant way and whether it is possible to have a non-trivial quantum contributions to the 
Energy-Momentum Tensor of Vacuum. 

We have addressed the first and in part the second issue on the basis of effective action method. 
In both cases the output of our investigation was perfectly consistent with the general expectations 
based on the known structure of renormalization in curved space-time and conservation law for the 
EMT. The calculations in the theory with SSB have shown some new term which was unnoticed 
until now. At the same time, after performing explicit calculations of this term we have found that 
the final results, Eqs. (|99p and (jlOOp . have the usual form and that the quantum effects always 
lead only to the renormalization of the inverse Newton constant and cosmological constant in Eqs. 
(f50l) and in the relations such as (ioTj) and (|52|) . 

It is important to note that our calculations can not be interpreted as a no-go theorem for a 
non-trivial quantum contributions to the low-energy sector of the gravitational action. As it was 
previously explained in [13], the chance to meet such corrections exists, but this can be verified 
only in the framework of some qualitatively new mathematical tool which should not be based 
on the perturbative expansion in curvatures. The linear in curvature approximation which was 
adopted here does not provide any information about these type of corrections. However, it was 
definitely worthwhile to check that the standard considerations really work for the non-trivial 
physical situations such as gravity combined with SSB. 
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Appendix A. Local momentum representation 



The calculations in Sect. 6 were done using Riemann normal coordinates 



see, e.g., 



4l[ for 



40J). In this Ap- 



introduction) and the local momentum representation technique (see, e.g., [39j, 
pendix we present some necessary elements of these tools and also derive the operator g^n — V^Vj, 
because it is related to relatively trivial calculations. 

The normal coordinates expansion performs around one special point P (that is related methods 
work well only for deriving local quantities), where metric is supposed to be flat Minkowski one. 
However, the derivatives of the metric, starting from the second one, are of course non-zero. The 
expression for the metric is 

v = V' ~ \ R m»p y a y 13 + ••■ ■ ( 101 ) 

Here and below all components of the curvature tensor correspond to the point P, also in (jlOlh we 
have omitted the higher order terms in curvature tensors and their derivatives. Furthermore y a 
represent deviation from the point P, such that all partial derivatives below are taken with respect 
to y a . It is fairly easy to derive, using (jlOip . the following expansions: 

= -~#V)^ r - ( 102 ) 

Then for the two covariant derivatives acting on scalar we obtain 

V^V U = 8^8 U + ^R X (fMu) T y T d x . (103) 

Making contraction with 

9" v = f + ^ M A^!/V, (104) 
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we get 



□ = <TV M V„ = d 2 - | R X T y T d x + l -R» a v p y a y^d^d u , 
where d 2 = rj^d^dy. Finally, the first operator of our interest is 

V M V^ -g^U = d^8 v + ^R x ( ^) T y T d x - ^ v 8 2 + -rj^R&dx 
- l^uR p a a y a y p d p d a - \R mvP y a yPd 2 . 



(105) 



(106) 



The next operator we are interested in is the propagator of scalar field. Direct calculations 
using (|105p (see, e.g., [3) Qj) l e& d to the relevant expression for the propagator of the field of the 



mass m in the linear in curvature approximation, 



G{y) = J 



d A k 
(27T)* 



1 



- k 2 + 



l 



R 



(107) 



2 V s 6/(/fc 2 +m 2 ) 2 - 

where we already assumed Wick rotation to Euclidean space. The Eq. ()107|) was recently used in 
32| to derive the effective potential of scalar field in the momentum cut-off regularization. Due to 
the use of local momentum representation (I107p the result has covariant form, despite the naive 
application of the cut-off scheme is supposed to break down even Lorentz invariance. 



Appendix B. Effective equations of motion 



Here we present the effective equations of motion on the (^-approximation, when the effective 
action is reduced to (|6ip . By using Eqs. (|61|) and (|66p we obtain 



and 



8 2 V n 



86 



8, 



V\ — R-FTT 



''() 



<-'() 



(108) 



(109) 



Let us denote the curvature-independent and mass-independent part of classical potential as 



V = V{6) = -6 



(110) 



From the quantities Vo and V\ given by (|63p and (|64p one can easily get 

tfV'" -(V"-m 2 )V"'ln^ 



Qydiv 


1 r 


86 


" 32vr 2 [ 


Q2ydiv 


1 r 


OS 2 


" 32tt 2 l 


Qydiv 


_ 1 ( 


86 32vr 2 V 



n 2 V"" - {V'"f In % - V""(V" - m 2 ) In % 



1 



n 2 



i--\V'"\n — 



(111) 
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Furthermore, from Vren and Ken in (1621) we obtain 



Qyren 

d<t> 

Q2yren 

d(f> 2 
Qyren 



1 



(V" - m 2 )V" 



32tt 2 
1 

32 



In 



V" - m? 



1 

+ 2 



Ki\(y'") 2 + (y" -m 2 )v" 



In 



V" - m< 



+ 



+ 



{V'"f} 



1 



In 



V" - m z 



90 32vr 2 V" 6/ L V jt< 

Next, we calculate the on-shell expressions by replacing 

1 



+ 1 



(112) 



and Xv 2 = 6m 2 , in the form 



Qydiv 



dV{ 



div 



dcj) 

Q2ydiv 



32n 2 L 

1 

32^ 
1 



XvqQ — 2 Am vq In — ^ 



m^ 



o/ m/ 

Aft 2 -8Am 2 In— ^ 
m 2 



t>o 327r 

Similarly, the analogous renormalized on-shell expressions are 
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At this point we can derive the elements of equations of motion, 
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(113) 



(114) 

(115) 
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(117) 
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The last observation is that, as we have already mentioned in the main text, all these expressions 
must be treated as constants in the given approximation. 
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